Experimental evidence of quantum quenching of fluctuations in heavy ion colhsions 
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The first experimental results of a new quantum method for calculating nuclear temperature 
and density of fragmenting heavy ions is presented. This method is based on fluctuations in the 
event quadrupole momentum and fragment multiplicity distributions of light Fermions. The cal- 
culated temperatures are lower than those obtained with a similar classical method. Quenching of 
the normalized multiplicity distributions of light fermions due to Pauli blocking is also observed. 
These results indicate a need for a quantum treatment when dealing with statistical properties of 
fragmenting heavy-ions. 
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The experimental measurement of temperature and 
density has significant impact on the field of intermedi- 
ate energy heavy-ion reactions. Constraining the nuclear 
equation of state and investigations of a liquid-gas phase 
transition in nuclear matter are two areas of intense in- 
terest that rely heavily on such quantities. Of broader 
interest, the nucleus itself provides a laboratory to look 
at the thermal properties of a quantum system. It was 
recognized some time ago by W. Bauer that it would be 
necessary to correct for the Fermi motion in probes de- 
rived from the kinetic energy of emitted fragments [ij. 
Despite this, most of the commonly used temperature 
calculations are derived from classical statistical mechan- 
ics 0,13] • One exception is a recently proposed calculation 
based on fluctuations in the event quadrupole momentum 
and particle multiplicity distributions [j, A similar 
method has been used with trapped Fermi gases, which 
shows a quenching in the density fluctuations of the gas 
due to Pauli blocking 0-01 ■ An analogous signal should 
be observable in the fragmentation of heavy ions. In this 
letter, we will present the first experimental quantum 
temperatures calculated by the method presented in Ref. 
Qualitative comparisons will be made to predictions 
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of the Constrained Molecular Dynamics Model (CoMD), 
a microscopic model that includes Fermionic statistics, 
for a similar fragmenting system [lo| . The experimen- 
tal effect of Pauli blocking on the particle multiplicity 
distributions of light fermions will also be presented. 

The experimental data was acquired at the Cyclotron 
Institute of Texas A&M University. A beam of 45 
MeV/A -^^S ions was accelerated by the K500 supercon- 
ducting cyclotron and then impinged upon a highly en- 
riched ^^^Sn target. Charged particles from the reaction 
were measured by the FAUST detector array. FAUST is 
a segmented array of 68 Si - CsI(Tl) telescopes with cov- 
erage from 1.6 to 44.8 degrees with respect to the beam 
axis. The array is structured into 5 square rings with 
the longitudinal distance from the silicon detector to the 
target kept constant within each ring. This arrangement 
gives FAUST solid angle coverage of 90% from 2.3 to 33.6 
degrees, 71% from 1.6 to 2.3 degrees, and 25% from 33.6 
to 44.8 degrees. A detailed description of FAUST may be 
found in Ref. [ill, [l2!| . It should be noted that the exper- 
imental setup was not designed to be sensitive to photons 
or neutrons; only charged particles were detected. 

Particle identification was performed using the AE-E 
technique. 



For details see Ref. [12 1 



Such identifica- 
tion requires the fragment to have sufficient energy to 
pass through the silicon detector before implanting in 
the CsI(Tl). Because of the thickness of the silicon, the 
charged particles collected are, by design, heavily biased 
towards fragments produced from projectile-like sources 
(PLS) in peripheral and mid-peripheral collisions. Par- 
ticles were identified isotopcially (in both atomic charge 
and mass number) for particles with a charge up to at 
least Z — 8 throughout the array and Z = 13 for select 
detectors. Particles for which the mass could not be de- 
termined were identified in charge only up to the charge 
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of the incident beam iZ = 16). 

PLS were reconstructed on an event-by-event basis. 
Events were selected based on two constraints. First, 
the sum of the charge of detected particles was required 
to be the charge of the beam (Z=16). This condition has 
previously been shown to select for events with properties 
indicative of a projectile-like source [H, 14 1. Second, the 



sum of the mass of the detected particles was required to 
be the mass of the beam (A=32). This also implies that 
all fragments in the event are mass identified and that 
the minimum charged particle multiplicity per event is 
two. By fixing the source to a single charge and mass, 
effects due to differences in the size or isospin asymmetry 
of the source should be significantly reduced. 

Event excitation was calculated by calorimetry. In or- 
der to minimize dynamical contributions to the kinetic 
energy 1^ 1^, only the components of the momentum 
transverse to the PLS-TLS (target-like source) separation 
axis were used in the calculation: 
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where K is the kinetic energy, and Py are the transverse 
components of a fragment's momentum in the frame of 
the PLS, and m is the fragment mass. The event excita- 
tion energy was calculated: 



AM,) - AMpLs 
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where K.^ is the kinetic energy of a fragment in the frame 
of the PLS, AMi is the mass excess of a given fragment 
and AMpLs is the mass excess of the PLS. The sum- 
mation is over all detected charged particles in an event. 
Events were sorted into 12 bins in event excitation energy, 
each bin has a width of 0.5 MeV. It should be noted that 
this estimate represents an apparent excitation since the 
experimental setup was not sensitive to emitted neutrons. 

The temperature calculation is based on the fluctu- 
ations in the event quadrupole moment. Details of this 
method can be found in reference 4] . Similar to the exci- 
tation energy calculation discussed above, any dynamical 
contribution to the fragment kinetic energy would largely 
affect the longitudinal component. The quadrupole mo- 
mentum is therefore defined as Qxy —< Px ~ Py > where 
Px and Py are the same transverse components of a given 
fragment's momentum used in the event excitation cal- 
culation. The variance of the quadrupole is then given 
by: 
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where f{p) is the momentum distribution of the particles. 
In earlier work, a classical Maxwell-Boltzmann distribu- 
tion was assumed for f{p) giving the variance as: 
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FIG. 1. Normalized multiplicity versus event excitation en- 
ergy is plotted for protons, tritons, and ^He. Notice that in 
all three cases the normalized multiplicity shows quenching of 
the multiplicity fluctuations compared to the classical limit 
(((AA'^)^) /N = 1) denoted by a dashed line. Statistical error 
is plotted for the normalized multiplicity and the excitation 
bins are 0.5 MeV in width. 



where Td is the classical temperature, m is the mass of 
the fragment, and N is the average multiplicity per event 
of the given fragment type [4, 15, 16]. This classical re- 
sult will be useful for later comparison to the quantum 
method. 

For the quantum calculation, the authors of Ref. Q 
constrained themselves to A=l fermions so we will do 
the same here. If a Fermi-Dirac distribution is used in 
Eq. 3, and £/ = e/o(p/po)^^/^^ = 36(p/po)^^/^^ MeV is 
the Fermi energy of nuclear matter, using an expansion 
on T/sf we obtain: 
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This gives us the quadrupole fluctuations in terms of tem- 
perature and density through e/ 0, S|- The value of 
e/o =36 MeV was chosen to be consistent with Ref. [3]. 

Following a similar derivation, the fluctuations of a 
fermion multiplicity distribution can be written as Q : 
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By substituting Eq. 6 into Eq. 5 the Fermi energy can be 
expressed in terms of quadrupole and multiplicity fluctu- 
ations. It is then simple to calculate the quantum tem- 
perature from either Eq. 4 or Eq. 5 It should be 
noted that the above derivation assumes that T/ej < 1. 
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FIG. 2. (Top) The proton temperature is plotted as a function 
of event excitation energy for both the classical and quantum 
calculations. The quantum method shows a significant re- 
duction in temperature over that of the classical limit for all 
excitation bins. (Bottom) The proton density as calculated 
through the quantum method is plotted versus event excita- 
tion energy. In both plots, statistical error is shown and the 
event excitation bins match that of Fig. [1] 



Fig. [T] shows the measured normalized multiplic- 
ity fluctuations (((AiV)^) /TV) for protons, tritons, and 
•^He particles as a function of excitation energy per nu- 
cleon. For a classical gas, the multiplicity of a specific 
cluster size will follow Poissonian statistics for which 
((A7V)2) /TV = 1 [7]. If, instead, we have a Fermi gas, 
the Pauli exclusion principal prevents multiple fermions 
of the same type from occupying the same state. This has 
the effect of blocking certain fragment emission channels 
leading to a reduction or quenching of the multiplicity 
fluctuations. Such quenching is observed in Fig. [T] for 
all three particle types and agrees qualitatively with the 
predictions of the CoMD model jHj as well as the results 
from experiments with trapped Fermi gasses 0, 01 . The 
increased quenching of the protons over that of the tri- 
tons, and '^He particles is believed to be due to the com- 
posite nature of the A=3 particles. Moreover, A=3 parti- 
cles have a smaller mean multiplicity compared to that of 
protons, and therefore A=3 particles have a smaller prob- 
ability for interacting with particles of their own type. 

Since the quantum derivation is specific to A=l 
fermions, only results based on the proton fluctuations 
will be presented from this point forward. A compari- 
son of the classical vs quantum temperature is plotted in 
the top panel Fig. [5] as a function of excitation energy 
per nucleon. There is a signiflcant increase in tempera- 



FIG. 3. Temperature versus density is plotted for both the 
classical and quantum temperature methods. A direct rela- 
tion can be seen with an increase in density corresponding 
to an increase in the difference in calculated temperature be- 
tween the two methods. 



ture in the classical over the quantum calculation with 
the difference exceeding 2 MeV in the highest excitation 
bins. The increased value of the classical temperature 
was predicted in references [H, 0] and is caused by clas- 
sical calculations failing to account for the portion of a 
fragment's kinetic energy due to the Fermi motion of the 
nucleons prior to emission. 

The bottom panel in Fig. [2] shows the proton density 
vs event excitation energy. The density was calculated 
assuming a value of po — 0.15 fni^^: 



(7) 



using the fermi energy calculated from Eq. 5 and Eq. 6. 
The range of densities, from approximately 0.01 to 0.03 
fm~'^ agree well with recent calculations by Wada, et al 

The density calculated is the local density felt by pro- 
tons emitted from the PLS. In a liquid-gas type phase 
transition, the protons represent the gas or low density 
region, which explains why the density is increasing with 
excitation energy. At low excitation it is reasonable to 
expect that most fragments are ejected through a sequen- 
tial surface type emission, i.e. at very low densities on 
the surface. As we move to higher excitation, more nu- 
cleons get excited from the liquid to the gas leading to 
a higher local density felt by the protons. Similarly, we 
would expect the density of the liquid to decrease with 
increasing event excitation. 
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FIG. 4. Energy density versus temperature is shown. The 
downturn in the data points around 2.5 MeV could indicate 
a region of rapid variation in the energy density indicative 
of a phase transition as predicted in Ref. 0]. However, the 
unavailability of lower temperature points makes this incon- 
clusive in this data set. Statistical error is shown. 



This evolution from surface emission to volume emis- 
sion may also explain features in both the normalized 
multiplicity and temperature plots. The effect of the 
Pauli blocking responsible for the quenching observed in 
Fig. [1] should be reduced with increasing event excita- 
tion since this will reduce the probability of two particles 
being at the same momentum. However, increases in den- 
sity increase the probability of interaction between parti- 
cles; thus these two effects are competing processes. This 
could explain why the normalized proton multiplicity in- 
creases only slightly from 2.5 to 7 MeV in event excita- 
tion. It is possible that a similar argument could be made 
for the slow onset of quenching in the case of the A=3 
particles. With increasing event excitation, the probabil- 
ity of emitting composite particles increases, leading to 
a greater chance of such particles interacting, and there- 
fore a greater effect due to Pauli blocking. However, the 
extension to such cases is not straightforward and will 
be the subject of future investigations. The effect of the 
density on the temperature is more easily seen by plot- 
ting these quantities against each other as shown in Fig. 
[3l The increase in the difference between the two temper- 
ature calculations with increasing density is easily seen, 
and is attributed to the increase in Pauli blocking with 
increasing density. 

Variation in the density could also explain the lack of a 
noticeable plateau or flattening in the event excitation vs 
temperature often referred to as the caloric curve. A flat- 



tening of the caloric curve is a feature of a system passing 
through a phase transition. This feature may be obscured 
for systems whose density is not constant over a range of 
temperature. For such as system, it may be instructive to 
compare the energy density versus the temperature. Fig. 
Slshows the energy density e =< E*y/A > p as a function 
of temperature. A similar result was obtained in Ref. Q , 
suggesting a liquid-gas phase transition near 3 MeV due 
to a rapid variation of the energy density at that tem- 
perature. Fig. m agrees qualitatively with the data from 
Ref. [1]; however, the rapid variation predicted falls at 
the lower limit of the experimental data. We stress that 
the system discussed in Ref. [4] was much larger (A=80) 
compared to the one selected in this work (A=32) and 
therefore only a qualitative comparison is possible. 

In conclusion, the first results of a new quantum 
method for the calculation of nuclear temperature and 
density have been presented. The temperatures calcu- 
lated show a substantial reduction when compared to 
a similarly derived classical method. This reduction in 
temperature is in agreement with previous theoretical 
predictions showing that classically derived methods over 
predict nuclear temperature by failing to treat the Fermi 
motion of nucleons. Quenching of the normalized multi- 
plicity due to Pauli blocking has been observed for three 
species of light fermions. Together these results highlight 
the need for a quantum treatment when dealing with sta- 
tistical properties of fragmenting heavy ions. 
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